Rigorous conditions for the existence of bound states at the 
threshold in the two-particle case.* 

Dmitry K. Gridne^ and Martin E. Garcia 
Theoretische Physik, Universitat Kassel, 
Heinrich-Plett-Str. 40, Kassel, Germany 

Abstract 

In the framework of non-relativistic quantum mechanics and with the help of the Greens functions 
formalism we study the behavior of weakly bound states as they approach the continuum threshold. 
Through estimating the Green's function for positive potentials we derive rigorously the upper 
bound on the wave function, which helps to control its falloff. In particular, we prove that for 
potentials whose repulsive part decays slower than 1/r 2 the bound states approaching the threshold 
do not spread and eventually become bound states at the threshold. This means that such systems 
never reach supersizes, which would extend far beyond the effective range of attraction. The 
method presented here is applicable in the many-body case. 
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I. INTRODUCTION 



In many problems of quantum mechanics it is important to know what happens to the 
wave function of a system as the bound state approaches the dissociation (decay) threshold. 
In particular, how does the size of the system in the ground state change as the system 
becomes loosely bound. Among multiple examples of loosely bound systems in physics one 
could mention negative atomic and molecular ions (1], Efimov states [2] and halo nuclei jjj, 

For a bound state of the system as it approaches the threshold there could be two possi- 
bilities. The first one is that the probability distribution given by this bound state spreads, 
meaning that the probability to find all particles together in the fixed bounded region of 
space goes to zero (the size of the system goes to infinity). One observes such spreading 
in helium dimer 4J or in halo nuclei like 6 He or n Li, which are so loosely bound that two 

neutrons are about to leave the system and form dilute nuclear matter around the core 

fl 

nucleus ( 4 He and 9 Li respectively) [3]. The second possibility is that the bound state does 
not spread and in this case it eventually becomes a bound state at the threshold (the size of 
the system remains finite). This phenomenon is called the eigenvalue absorption. This is 
the case for doubly negative ions and proton halos [l] , [3] - 

Recall, that for two particles interacting through spherically symmetric potentials with fi- 
nite range S— states always spread, while all states with non-zero angular momentum become 
bound p]. Incidentally, it is natural to conjecture that the ground state of a multi-particle 
system with pair interactions of finite range (Vij(x) = for \x\ > R) cannot be bound at the 
threshold given that the particles are either bosons or distinguishable. For fermions with 
short-range interactions it is hard to say from general principles whether the ground state 
would spread or not. The physical approach in this case is to use some kind of shell model 
and to figure out if there is a centrifugal barrier, which prevents the wave function from 
spreading. 

On the other hand, there are potentials, for which bound states do not spread at all and 
when approaching the continuum they give rise to bound states exactly at the threshold [3] , 
js]. In particular, the physically important case of repulsive Coulomb tail case belongs to 
this type (see the discussion in |9[). Let us illustrate this situation by a simple example. 
Consider the square well potential plus a repulsive Coulomb tail, as in Fig. 1 (left), and 
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imagine the ground state in this potential as it approaches the threshold. The probability 
distribution for this state would remain a confined wave packet regardless of how small is 
the binding energy. For zero binding energy there would be a bound state, which would 
have a falloff of the type ~ exp(— \/r). On the contrary, if we cut off the positive tail at 
some arbitrary distance R c , the ground state approaching the threshold would eventually 
spread when the binding energy is sufficiently small, i.e. the probability to find the particle 
in some bounded region of space goes to zero with the binding energy. The state "tunnels" 
through the barrier. Note that this change in the behavior does not depend on the value of 
R c , which could be as large as pleased, so this effect is solely due to the repulsive Coulomb 
tail. This requires a special care in numerical calculations of loosely bound systems, because 
often in the calculations the potential becomes effectively cut off. 

A rigorous proof of the eigenvalue absorption in the case of a general short-range potential 
plus a repulsive Coulomb tail was given in Q . Unfortunately, the approach presented in , 
based on the Green's function expansion, is aimed specifically at the Coulomb long-range 
part and does not allow generalizations, for example, to potentials having long-range parts 
of the form r -1 + Ar~ 2 , which may arise in multipole expansions. Here, we shall consider 
a more general two-body case. The analysis in 7| illuminates the possibilities for various 
radially dependent long-range parts. However, the arguments given there are not rigorous. 
Finally, we should mention, that the phenomenon of the ground state absorption was proved 
rigorously in 10] for a three-body system with pure Coulomb interactions and the infinitely 
massive core. This makes one conjecture that at the point of critical charge negative ions 
have bound states at the threshold, see the discussion in In a forthcoming article we 
shall give the rigorous proof of this conjecture for distinguishable particles and bosons 11], 
see also the discussion concerning many-body systems in the last section. 

The paper is organized as follows. In Sec. [Til we set the criterium for the eigenvalue 
absorption. In Sec. II I II we derive useful upper bounds for the Green's function. These 
bounds can also be used to control numerical solutions. In Sec. UV]we prove our main result 
saying that potentials decaying slower as 1/r 2 give rise to bound states at the threshold. 
The last section presents conclusions and a short discussion concerning many-body systems. 
Finally, the Appendix contains technical details necessary for the proof in Sec. [TTJ 
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FIG. 1: Sketched behavior of the ground state wave functions approaching the threshold in two 
potentials. Left: the potential has a positive Coulomb tail. Right: the same potential is cut off at 
some distance R c . 

II. BOUND STATES NEAR THRESHOLD 



In nature we can make the bound state of the system approach the threshold by chang- 
ing the number and the type of particles. In theory we reproduce this behavior changing 
continuously some parameters in the system. For example, in the case of ions diminishing 
the atomic charge Z down to the critical value Z = Z cr makes the ground state approach 
the threshold l]. Here, for parameter, changing which we force the states approach the 
threshold, we take the coupling constant of the interaction. 

In our analysis we shall consider the Hamiltonian of two particles H = H + XW, where 
Ho = p 2 is the free Hamiltonian (we use the units where h = 1 and m = 1/2), W is the 
interaction and A is a coupling constant. By decreasing A we can lift any bound state to 
the continuum. For convenience we shall consider only W G L 2 (R 3 ) + L°°(M 3 ) 12] . This is 
a large class of interactions, which allows singularities not worse than r~ a for a < 3/2. In 
this case H is self-adjoint on the domain D(H ) = D(—A) fl2j |. However, one could extend 
our results to potentials having singularities of the type r~ a for a < 2. 

Let us assume that there is a bound state having the energy E(X) < for some value 
of the coupling constant A . E(\) increases monotonously when A decreases and eventually 
E(\) becomes zero for A = A cr , where A cr is called the critical coupling constant p|. In 
the following we show, using a simple example, how the wave function spreads in the case 
of exponentially decaying potentials. In this simple case the analyticity of the energy as a 
function of A helps us to establish an upper bound on the wave function 
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Theorem 1. If there exist A, a > such that \W\ < Ae~ a ^ and at X = X cr there is no 
zero- energy bound state then the following upper bound holds for the bound state ip having 
the energy E(X) in the neighborhood of E(X cr ) = 

M<— , 1 

r 

where C > is some constant. 

Proof, ij) satisfies the integral equation if) = —{H — E]~ l Wip, which could be rewritten as 



J 4n\x — y\ 

Using \W\ < Ae~ a ^ and applying the Schwarz inequality to Eq. (j2J) gives us 



\ip\ < A(^||^||^) 1/2 



j{ e -o-\y\ e -2y/\E\\x-y\ 

d y - — — 2 



1/2 



<\mw\\i)) i/2 C ' e ~ lElr (3) 



F - y\ 

\\'}wvc (" is some constant and r = \x\. On the other hand. i/ccal] (31 that at A = X cr 
the energy E(X) is analytic and can be expanded into convergent power series E(X) = 
Y^=2 a k(^ — Kr) k , where a<i < (a\ = because by condition there is no zero-energy bound 
state at A = A cr ). Applying the Hellman-Feynman theorem (?p\Ho\ip) = E(X) — XdE/dX 
gives us (iP\Hq\iP) /IE] 1 / 2 = 0(1). Because \W\ < Ae~ a ^ there must exist such constant L 
that \W\ < L(2r)~ 2 and thus by the uncertainty principle Q, [l3] (0||W||0) < L (<f>\H \<f>) 
for any G D(H ). Thus (V'HWHV') /\E\ 1/2 = 0(1), which together with Eq. j3]) proves the 
statement. □ 

As one can easily see, the function on the right-hand side of Eq. (CQ) dominates the wave 
function and spreads as E — >• maintaining the constant norm independent of E. The 
probability to find the particle in some fixed region of space goes to zero. 

Now, let us consider potentials with positive tails. Throughout the paper we shall assume 
for such potentials that W(x) > for \x\ > Rq, i.e., that outside some sphere the positive 
part dominates. Below we present a simple criterion, which tells us when bound states 
become bound states at the threshold. From the discussion above and from Theorem [1] it 
is clear that the strategy could be proving that as A \ X cr bound states remain confined 
in some region of space, i. e. do not spread. One way to achieve this is to show that bound 
states are dominated by some fixed function g 6 L 2 . 



Theorem 2. Let X n be a sequence of coupling constants and X n \ X cr . The following is 
true (a) if for each X n there is a bound state <fr n with the energy E n such that \<fi n \ < g, 
where g G L 2 and E n — > then there exists a normalized bound state at the threshold, that is 
4>o G D(H) and H(X cr )(f)Q = 0. (b) if for each X n there are m orthogonal bound states (j)^ 
with the energies E^ such that \<pn \ < g, where g G L 2 and lim^oo En = then there 
exist m orthonormal bound states at the threshold. 

The proof of this theorem, which follows the method in jl^ j. is given in the Appendix. 
Intuitively, this proposition is obvious. If the states do not spread they should finally 
form some bound state at the threshold. A similar theorem concerning many-body systems 



appeared in Zhislin and Zhizhenkova 15]. There the authors proved that if a minimizing 
sequence for the energy functionals does not spread, then there exists a minimizer in L 2 . 
Their result could be explained from practitioners point of view. Imagine that the system 
has no bound states with negative energy. Then, if the function minimizing the energy 
functional, does not go to zero as the number of basis functions increases, then there must 
exist a zero-energy bound state. 

Now let us see how the criterion in Theorem 2 works. First, we separate positive and neg- 
ative parts of the potential W = W + — W-, where W + = max(0, W) and W_ = max(0, — W) 
and W± > 0. The equation for the bound states reads H(X)<p = —k 2 <p, where k — > as 
A — > X cr . This can be rewritten as 

(H + k 2 + XW+)(j) = AW_0 (4) 

or equivalently 

= A(# + k 2 + XW^W^cf) (5) 
The operator (H + k 2 + XW + )~ l is an integral operator, whose kernel is positive and real 



171 ] . Thus we can rewrite Eq. (jSJ) as 

101 < 2A cr (# + k 2 + AW+)~ 1 W'_|0|, (6) 

where because A \ X cr we have taken A < 2A cr without loss of generality. If we show 
that the right-hand side of Eq. ([6]) is bounded by some fixed square integrable function, 
then according to Theorem [2] we would have bound states at the threshold. The operator 
(Hq + k 2 + XW+)" 1 is an integral operator, and its kernel is the Green's function having 
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two arguments. Because the function vanishes outside some sphere, the behavior of 

|0| at infinity is determined by the asymptotic of the Green's function when the integration 
argument is fixed within the sphere. Thus to find the asymptotic we need to derive upper 
and lower bounds on the Green's function. 



III. BOUNDS ON THE GREEN'S FUNCTIONS 
A. Potential tails decaying as 1/r 2 
1. Upper Bound 

We introduce the function which would play the role of potential's tail: 

f if r < R 

V (A,R 0] x) = { (7) 
[ Ar~ 2 if r > #0, 

We are interested in the kernel of the integral operator [Hq + k 2 + rj(A, Rq, for k 
real, which we denote as Gk{A, R ; x, y). Note that the kernel of such operator is a positive 



function, continuous away from x = y 17j]. Our aim in this section is to find the upper 



bound on Gk(A, Rq; x, y). For that we need the following lemma 

Lemma 1. Let Gi^ 2 (x,y) denote the integral kernels of [H + k 2 + Vi^] -1 and suppose 
V 1 (x)<V 2 (x). ThenG 2 (x,y)<G 1 (x,y). 

Proof. Through the integral representation we get 

G 2 (x,y) = Gfay) - [ dy / G 1 (x,y')[V 2 (y') - V x { y ')}G 2 {y' ,y) (8) 



Eq. (jSj) is the kernel representation of the equation (A + B)^ 1 = A^ 1 — A^ 1 B(A + B)^ 1 , 
where A = H + k 2 + V\ and B — V 2 — V%. Now, because Gi, 2 (x,y) are positive [17] and 
the potential difference V 2 — V\ is non-negative, the integral in Eq. (jSj) must be non-negative 
and the statement is proved. □ 

There is another elucidating and more direct way to prove Lemma [U The proposition of 
the Lemma follows from the Laplace transform and the Trotter product formula, which lie 
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at the heart of path-integrals [1^ 

(i?o + k 2 + V i y 1 = I e- kH e^ Ho+v ^dt (9) 



e 



o 

-t(H +V!) _ „_ i- l„-tH /m„-tVi/m 



,-lim [e-tHo/mg-m/mj'" (10) 



Because e~ tH ° in Eq. fflUj) has a positive kernel, namely (47it)^ 3 ^ 2 e~^ x ~ y ^ ^ u , the kernel of the 
operator on the left in Eq. (TTOT) becomes smaller when Vi is replaced by V^- 

The idea behind the upper bound on the Green's function is rather simple. Suppose 
we have found such a function F(A, Rq, x), independent of k, so that G k (A, R ; x, 0) < 
F(A, Rq; x) holds for all k and x. First, we shall derive the upper bound in terms of the 
function F, and then we shall determine F explicitly. 

Let us fix the functions A(s), Rq(s) so that the following inequality holds 

7 ] {A,R ;x)>r ] {A,R ;x-s), (11) 

where s is some fixed three-dimensional vector. From simple geometric arguments it follows 
that Eq. (flT!) would be satisfied if A(s), Ro(s) satisfy the inequalities 

Ms) >R + \s\ (12) 
A(s) < A^-^ (13) 

Let us mention that the closer A is to A the better is the asymptotic behavior of the bound, 
hence it is reasonable to take Rq large. 

Translating the arguments one finds that Gk{A, R ; x — s,y — s) is the integral kernel of 
the operator [Hq + k 2 + rj(A, Rq; x — s)] _1 . Now, using Eq. ( TlTT) and Lemma [1] we obtain the 
upper bound 

G k (A, R ; x, y) < G k (A, R ;x-s,y- s) (14) 
Eq. ffT4|) is valid for all s, so we can put s = y, which gives us 

G k (A, R ; x, y) < F(A(y), R^{y); x - y) (15) 

It remains to find F, which is the upper bound on G k (A, Rq; x, 0). This is easy because 
G k {A, Rq; x, 0) is spherically symmetric in x. From now on for simplicity of notation we 
shall drop A, R in the arguments, writing, for example, F(x) instead of F(A, R Q ; x). First, 
we shall give a formal solution, then we shall prove that this solution is indeed correct. By 
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Lemma [T] 0) < Gk>(x, 0) if k! < k, so we can take F(x) = lim^ G k (x, 0). Because 
Gk(x, 0) is continuous away from x = [171] and the functions Gk(x, 0) increase monotonically 
when k — > the pointwise limit makes sense. By Lemma [T] G k (x, 0) < Crjj. (x, 0), where 
G k °\x,y) = (4ir\x — y|) _1 exp (— k\x — y\) is the free propagator, i.e. the integral kernel of 
the operator [Ho + k 2 ]^ 1 . This means F(x) is bounded away from x = and _F(x) < (47rr) -1 . 
Because G k (x,y), formally satisfies the equation [H + k 2 + rj]Gk{x, y) = 5(x — y) one expects 
that F(x) satisfies the equation 

[H + V }F = 5(x) (16) 
To find the solution of Eq. (FlBl) we set 



1 1 + br if r < R 



x 



47rr cr~ a if r > #0, 



(17) 



where a is the positive root of the equation a(a + 1) = A and the constants b, c are fixed 
requiring, as usual, that F and its derivative are continuous at r = R . This gives us 

1 f 1 - i?n Va + 1)-V ifr<i?o 
F(AjRo;x)= X J ° \ ^ (18) 

4?rr I + a)-V- a ifr>.Ro 

One can check that F(x) defined by Eq. ([TBI indeed satisfies Eq. fTTBT) . 
For completeness we give the accurate proof, which justifies Eq. ( FlBl . 

Lemma 2. F(x) defined as F(x) = lim^ G k (x, 0) equals a.e. the expression given by 
Eq. m). 



Proof. The integral equation for the resolvent reads |]J 

G k {x,y) = Gf (*,y) - [ dy'G k °\x iy ') V (y')G k (y\y) (19) 



Substituting the expression for G k (x, y) and setting !/=0we obtain 

Gk (x,0) = ^—-- dy'- jrvWGM, 0) (20) 

47rr 47T J \x — y \ 

Applying lim^o to both sides of Eq. (120]) gives us the integral equation 

= T~ J~ / ^'t^V^/) (21) 

47rr 47r J |x — y'| 

By simple substitution and calculating the integrals one can check that F given by Eq. (fl8|) 

indeed solves the integral equation Eq. (121"]) . It remains to prove that no other solution 



exists. Suppose there are two solutions and denote their difference Z = F\ — F2. Then Z 
satisfies the integral equation 

Z W = / dy'^^{y')Z{y') (22) 

We need to show that Z = a.e. Let Z = if^Z, then Z E L 2 by the dominated convergence 
theorem and \\Z\\ 7^ 0, because otherwise from Eq. (1221) it follows Z = and we are done. 
From Eq. ([22]) we obtain 

k — 2/I 

But Eq. ( 1231) cannot hold because \x — y\~ l is the kernel of the strictly positive operator. 
This means Eq. (1221) holds only if Z = 0. □ 

Now let us formulate the bound in the form required in Sec. IIVI 

Corollary 1. Let A > 3/4, then there exist C > 0, 5 > and R such that for \y\ < R Q and 
\x\ > R the inequality holds Gk(A, Rq; x,y) < C\x\~ 3 / 2 ~ s . 

Proof. For \y\ < Rq we can fix the values Rq and A independent of y. Both inequalities 
Eq. (fT2lh (TI3I) would be satisfied if .Ro > 2i? and A < AR 2 {FLo + R )~ 2 . When .Ro becomes 
large A gets close to A, so we can fix the values of Rq and A to ensure that the following 
inequality holds 3/4 < A < A. If we set R = Rq + Rq, then for \y\ < Rq and \x\ > R we 
have \x — y\ > Rq and from Eq. (TT5|) . (fl8l) we get 

G fc (A, i?o; x, y) < (4tt|x|)- 1 ^(1 + ^{x - y\~~ a < Clx)" 1 " 5 , (24) 

where a is the positive root of the equation a{a + 1) = A, a > 1/2. □ 

J?. Lower Bound 

Here we shall briefly discuss how the same method can be applied to construction of 
lower bounds. We would need this in Section IIVI where we show that the ground state of 
potentials decaying faster than (3/4)r~ 2 spreads near the point of critical binding. For that 
we need the following type of potential 

, V if r < Ro 

Z(A,Rq,;x)={ (25) 
Ar~ 2 if r > R , 
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and we need the upper bound for Green's function of the operator Ek(A, Ro) = [Ho+^+k 2 ]^ 1 , 
which has the integral kernel Sk(A, Ro] x,y). We shall derive the lower bound in terms of 
the function fk(A, R ; r) = E k (A, Rq; x, 0), which falls off at infinity and solves the following 
equation 

[H + £ + k 2 ]f k = 6(x) (26) 

fk depends only on r = \x\ (because the potential is spherically symmetric) and is a con- 
tinuous function away from r = 0. By definition of fk we have < fk < l/(47rr). Setting 
fk = ( / lirr)~ 1 fk from Eq. fT26|) we obtain the equation on fk 

~fk+ tfk + k 2 h = (27) 

with the boundary conditions jfc(0) = 1 and /fc(oo) = 0. The function < fk(r) < 1 comes 
out as a solution of a simple radial equation and thus can be easily calculated. As usual, one 
calculates the solutions fk{r < Rq) and fk{r > Ro) and determines the constants so that fk 
and its derivative are continuous at Rq- The following Lemma is useful for the lower bound. 

Lemma 3. For r > R there exists Co independent of k such that 

fk(r) > C e- kr r- a (28) 

Proof. According to Eq. ( 1271) on the interval [Rq, oo] the function fk satisfies the equation 

-f" + Ar- 2 fk + k 2 fk = (29) 

Let us set fk{ r ) = gk(r)e~ kr r~ a . Then for gk(r) on the interval [Rq, oo] the equation becomes 

-gfi + 2 {ar' 1 + k)g' k = 2ak r - 1 g k (30) 

Because fk is positive gk should be also positive. Hence from Eq. fl30|) we get 

g'k' < 2 (ar- 1 + k) g' k (31) 

We want to show that g' k > 0. Indeed, if on the contrary g' k (y) < at some point y, then 
at this point due to Eq. (13"T|) g' k \y) < 0. Hence g' is a monotonically decreasing function for 
r > y. Thus from Eq. ( 13TT) we conclude that g' k < 2kg' k {y) for all r > y, i.e. the second 
derivative is less than a fixed negative value, which means that at some point gk becomes 
negative. Hence the assumption was false and g' k > holds. On the other hand, fk > 
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and as k — > the function f k monotonically increases in all points. Hence, there must exist 
Co > such that gk(Ro) > Co. Together with g' k > this means that g k stays above Co and 
Eq. (PSD holds. □ 

Now we follow the above procedure and define A, i? as satisfying the inequality 

C{A,Ro,;x-s) >£{A,Ro,;x) (32) 

By geometrical arguments A, R must satisfy 

Ro(s) > R + \s\ (33) 

Ms) > A g (34) 

Just as in the previous subsection through Eq. (|32|) we obtain the desired lower bound 

S fc (A,i2o;x,y) > E k (A,R ;x - y,0) = f k (A,R ; \x - y\) (35) 

Now suppose A < 3/4. Looking at Eq. (|33l)-(l34l one can see that we can fix Rq and A 
so that A < 3/4 and Eq. ( |33l) - (|34l) are valid. Then from Eq. (135j) and Lemma [3] it is clear 
that there exists a constant C > such that for |x| > 2i? and for \y\ < Rq the following 
inequality holds 

E k {A, R ; x, y) > f k (A, R - \x-y\)> Ce^lx^ 2 (36) 
We would need inequality Eq. (136!) in Sec. HVl 

B. Potential tails decaying as 1/r 

Here we would like to apply the results of the previous section to potentials with positive 
Coulomb tails. This helps to establish the decay properties of eigenfunctions lying at the 
threshold. We shall not present a detailed exposition, because everything is similar to the 
previous section. One can follow in the steps of the previous section and derive the bound in 
terms of the solution of the equation [Hq + rf\F = S(r), where rj' is the Coulomb tail. This 
however could not be expressed through elementary functions, so we shall make a couple of 
simplifying approximations. We shall consider the following potential tail. 

,0 if r < R 

C(a J R ;x)={ (37) 
{a 2 /A)^ 1 + (a/4)r~ 3 / 2 if r > Rq, 
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The repulsive Coulomb tail dominates in the potential of Eq. (1371) and one can choose 
the constants so that the actual Coulomb tail is greater than the function in Eq. fl37j) . Let 
G c k (a, R ; x, y) be the integral kernel of the operator [H + k 2 + ((a, R ; x)]" 1 , where c stands 
for Coulomb. The rest follows as above. 

Let us fix the functions ci(s),Rq(s) so that the following inequality holds 

((a,R ;x) > ((a,R ;x - s), (38) 

Again from geometric arguments it follows that Eq. (|38|) would be satisfied if a(s),Ro(s) 
satisfy the inequalities 

Ro{s) > R + \s\ (39) 

~a(s) < a (J^j-) (40) 
\R + \s\ J 

where the above conditions are obtained by directly applying Eq. (1381) to each positive term 
in the expression for ((a, R ; x) given by Eq. (1371) . 

Again let us define F c (x) = lim fc ^ Gl(x, 0), which makes F c (x) satisfy the equation 

[H + QF c = 6(x) (41) 

As one can easily check, the solution of Eq. (I4TT) is given by 

1 I 1 - __i__ r if r < R Q 

F c (A,R 0] r) = x<^ R^^m/a - ^ 

4vrr I e»v^ -oVr if r > d k j 

{ l+(a/2)v / % C U ' ~ 110 

Finally, the upper bound reads 

G c k (a, R ; x, y) < F c (a(y), R (y); x-y), (43) 

where a and Rq satisfy Eq. fl39|) - fj40|) . As in Corollary [1] from Eq. fj43|) we find that there 
exists such R > and C > that 

G c k (a, R ; x, y) < Ce~ a ^ if \y\ < R and \x\ > R (44) 

As we have mentioned for potentials with positive tails the asymptotic of the Green's function 
determines the fall-off behavior of bound state wave functions. Hence the bound state wave 
functions fall off at least as fast as e~ a ^. Calculating in the same way the lower bound one 
finds that this is the actual fall-off. 
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IV. MAIN RESULT 



Now we state the main result of this paper. 

Theorem 3. // there are such Ro and A > 3/4 that XW + > i](A, Rq; x) then at A = A cr all 
states that hit the threshold at A = A cr become zero energy bound states. 

Proof. Let us define Gk(A, R ; x, y) the positive integral kernel of the operator [H +k 2 +riY l . 
Then from Eq. and by Lemma [TJ we get the bound 

x) < 2\ cr [ dyG k (A,Ro;x,y)W-(y)\<l>\{y), (45) 
J\y\<Ro 

where we have used W^{y) = for \y\ < i? - Now we shall use the upper bounds on the 
Green's function Gk(A, R Q ; x,y) derived in Se. IHII For \x\ > R we can use Corollary [I] to 
obtain from Eq. fj4"5]) 



x) < 2\ cr C\x\-^ 2 - s / dyW-(y)\<l>\{y) < C^xl^ 2 - 6 = 9> (x), (46) 
Ay\<Ro 

where we have applied the Schwarz inequality and used W e L 2 + L°°. For \x\ < R we can 
use Gk{A, R ; x, y) < (47r) _1 |a; — y\~ x to obtain from Eq. fj4"5l) 



x) < 2A cr (4vr)- 1 / dy\x- y\- x W-{y)\<l>\(y) 



< G 



dy\x-y\ 2 W 2 (y) 
U\y\<R 



1/2 

= 9< (x), (47) 



Thus we get \(f>\(x) < g(x), where g{x) = g<(x) for |x| < R and g(x) = g>(x) for |x| > R. 
Because g(x) G L 2 Theorem [2] applies and the theorem is proved. □ 

One can construct simple examples, which show that the condition in Theorem [3] is best 
possible. The following theorem is also true. 

Theorem 4. // there exists Rq such that 

W(x) < (3/A)\x\- 2 forlxl^o (48) 

then the ground state cannot be bound at the threshold. 

Proof. For simplicity we shall assume that there exists a constant Vq such that W + < Vq. A 
proof by contradiction. Suppose that the ground state tp exists. Then the equation for the 
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bound state at the threshold can be written as [Ho + W + ]ipo = W-ifro. This is equivalent to 
the equation [Hq + W + + k 2 ]ipo = W-ip + k 2 tpo, which in turn can be transformed into the 
integral equation 

ij = [H + W+ + k^W.tPo + k 2 [H + W+ + /c 2 ]~Vo (49) 



The ground state wave function is always nonnegative ip > [12|. Because W_ > 0, 
ipo > in Eq. ( B9i) on the right-hand side we have a sum of two positive terms (the operator 
[H + W + + k 2 ]~ x has a positive integral kernel). Hence we must have 

|| [H + W+ + A; 2 ]- 1 ^.^!! < 1 (50) 

for all k. Because the positive part of W is bounded we have W + < £(A, R ; x), where ^ is 
defined by Eq. (|25l). 

From Eq. fl50l) and using Lemma [1] we conclude ||Hjt(^4, i?o)Wl-"0o|| — 1 for all k, where 
E k (A,R ) = [H + (, + k 2 }^ 1 . Our aim is to prove lim fe ^ ||Hfe(A, i2 )W_^o|| — 00 thus 
obtaining the desired contradiction. We shall use the lower bound on Sfc(/L, Rq; x, y) from 
Sec. IIII A 21 Let us fix Rq as in the last part of Sec. IIII A 21 Then using the bound Eq. (1361) 
we obtain for the square of the norm 



\\E k W^ \\ > / / dy 1 dy 2 W-(v 1 )MyiW-(m)My2) (51) 
Avi\<Rq Av2\<Ro 

x / dxE k (A,R Q ;x,y 1 )E k (A,R ;x,y 2 ) (52) 

J\x\>2Ro 

> M 2 C 2 [ dx\x\- 3 e- 2k W (53) 

J\x\>2R 

where M = J^ <RQ dyW_(y)ilJo(y) is some fixed constant. Note that M^O because that 

would mean ip = 0. It is clear that the right-hand side in Eq. (I53p becomes infinitely large 

as k — > and thus i^W-^oll — 1 cannot hold for small k. □ 



V. CONCLUSIONS 

We have proposed the method to derive lower and upper bounds on the Green's functions, 
which helps to determine the fall-off of bound states. Using these bounds we have proved 
that potentials, whose tails decay as {2ji/h 2 )V > (3/4)r~ 2 , where \x is the reduced mass, 
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FIG. 2: Typical stability diagram (sketch) for three Coulomb charges {—lj 91,92}, the shaded area 
representing stable systems. On the arcs of stability curve where either q\ > 1 or q 2 > 1 there are 
bound states at the threshold. 

absorb the eigenvalues, meaning that their bound states do not spread and become bound 
states at the threshold. We have also found that ground states in potentials, whose tails 
decay as (2fi/h 2 )V < (3/4)r~ 2 , always spread as they approach the continuum. 

These methods can be applied to the many-particle case, where it is still not known, 
for which pair interactions between decaying particles or clusters the bound state would 
become absorbed. The difficulty is that it is hard to control the asymptotic of the many- 
body bound state wave function. Using bounds on the Green's functions derived here one 



can demonstrate 



ll| that when there is a long-range Coulomb repulsion between decaying 



components (particles or clusters), then the bound state must get absorbed. 

Two types of behavior, namely spreading and eigenvalue absorption can be perfectly 
illustrated by a stability diagram for three Coulomb charges, see 18[. When masses are 
fixed the diagram has the form as in Fig. 2. If we consider, for example, the upper arc, 



which is the stability border, it has the following property 18j. Up to some non-zero value 
of the stability border is given by the equation q 2 = 1. Then at the point {q®, 1} the 



arc goes up. With the same method as here it can be proved [TJJ that if one approaches 
the stability border from the side where q 2 = 1 then the ground state spreads and there 
is no bound state at the threshold. On the contrary, for the points on the arc the ground 
state becomes absorbed, i. e. it does not spread and becomes a bound state at the threshold. 
The reason is that on the stability border, where q 2 = 1 the system decays into a neutral 
cluster and one charged particle, and for q 2 > 1 both cluster and the particle are charged 
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positively. The resulting Coulomb repulsion between these objects hinders the spreading of 
the wave function and the ground state becomes absorbed. Note that it has already been 



proved rigorously [10| , that in the case of an infinite core and two other masses being equal, 
the sharp point on the diagram, see Fig. 2, has a bound state at the threshold. Our method 
helps to extend this result to many particles. 

APPENDIX A: PROOF OF THEOREM [2] 

Proof. Let us prove part (a). We follow the argument from [14!]. Because \\(f) n \\ = 1 we can 
extract a weakly converging subsequence (for which we reserve the same index n) such that 
<fin —> <fio, where <fio £ L? . Because H is self-adjoint in order to prove that <fi§ £ D(H) and 
H(\ cr )(fi = it is enough to show that for every / £ D(H) we have (H(X cr )f\(fi ) = 0. The 
latter we obtain as follows 

(H(X cr )f\<fi ) = lim (H(X cr )f\<fi n ) = lim [(H(X n )f\<fi n ) - (X n - X cr ) (Wf\<fi n )} = (Al) 

lim (f\H(X n )<fi n ) = lim E n {f\<fi n ) = 0. (A2) 

n^oo n^oo 

The only thing that remains to show is that 0o 7^ 0. We shall prove this by contradiction 
assuming that <fi n 0. Let us introduce xr-, the characteristic function of the interval 
[0,.R] {i.e. Xr( x ) — 1 when \x\ £ [0, R] and Xr( x ) — otherwise). Because \<fi n \ < g and 
g £ L 2 we can fix R so that {(fi n \XR\ < fin) > 1/2. We would like to show that for <fi n — > 
the condition (<fi n \XR\<fin) > 1/2 cannot hold for large n. One way to do this is to use 
that J \ V(fi n \ 2 dx < const and apply the Rellich-Kondrashov lemma [l6f] giving XR'fin 
strongly. Or we can use the argument similar to the one in [l4]. Using the equation 
(H + l)(f> n = (E n + l)0 n - X n WK we get 0„ = (E n + 1)(H + 1)"V„ - A n (^ + l)~ l W<fi n . 
Substituting this into ((finlXRl&n) > 1/2 we obtain 

(E n + 1) ((fin\XR(Ho + 1)-Vn) - X n {<fi n \x R (H Q + l)" 1 ^) > 1/2 (A3) 

The operators xr(Hq + ^)~ 1 and xr(Ho + 1)~ 1 W have square integrable kernels and are there- 
fore compact. Acting on weakly convergent sequences they make them converge strongly 
and hence both terms on the left-hand side of Eq. (IA3I) go to zero. Thus Eq. fl A3j) cannot 
hold for large n, which proves (a). 

Part (b) easily follows if we prove that from <fi n — > 0o follows <fi n (fio in norm. Indeed, 
for each A n there are bound states (fin (i = 1, m) satisfying H(X n )(fi$ = E n ^(fin\ 
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Moreover \<p n j — So- and as X n \ A cr the energies go to zero En — > 0. In this case 
since \(j)n\ < g <E L 2 there are m bound states at the threshold, (j)n — > (fy£' . Because this 
convergence is in norm ^o^o / = holds. 

To prove that from <p n — > <fi follows 4> n — > 0o in norm let us define £ n = 0„ — 0o, then 
£ n — > and we would like to show that ||£ n || — * 0. A proof by contradiction. If not then 
there must exist a constant a > and a subsequence (for which we again reserve the same 
index n) such that ||£ n || 2 > Again because \(j) n \ < g and <po G L 2 we can fix R so that 
(Ux*l£n) > a/2. We have f n G D(H) and (if + A„W)£ n = E n (j> n + (A cr - A n )lT0 o - From 
this equation we easily get 

U = (H + l)" 1 ^ - A n (# + l) _1 ^n + £n(#0 + l)~Vn + (A cr - A n )(#o + l) -1 ^ (A4) 

Substituting one £ n from Eq. (1A4jl into (£ n , Xi?£n) and using that xr(-^o + I) -1 and xr(Hq + 
1)~ 1 W are compact and ( n ^ we obtain (£, n \XR\£,n) 0. This is a contradiction. □ 
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